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1 Introduction 

Diffraction experiments are important to determine the structure of a solid, even more so with 
the refined methods available today. Recent applications include aperiodic systems as well as 
systems with disorder. 

Kinematic diffraction can be described and understood in terms of Fourier analysis. The 
diffraction image is related to the Fourier transform of the autocorrelation or Patterson function 
of the scattering obstacle, e.g., the electron density. The situation is well understood for (ideal) 
crystals, which show a complete lattice of periods, even though the corresponding inverse 
problem does not have a unique solution in general. 

Beyond the periodic situation, firm results are sparse, and until recently, one did not even 
know which general distribution of scatterers would result in a pure point diffraction spectrum, 
i.e., in a diffraction image of pure Bragg peaks only. This review revolves around this ques- 
tion m, and summarises the present state of affairs, with special emphasis on contributions 
obtained during the time of the DFG focus program on quasicrystals. We concentrate on the 
few results here that have been established rigorously, and indicate where further research is 
needed. 

The article is organised as follows: In the first section, we state conditions (TheoremQ) 
under which a set of scatterers diffract, i.e., its diffraction spectrum consists of Bragg peaks 
only. Instead of looking at very general set of scatterers, we specialise in those sets that 
model physical structures (as weighted Dirac combs). Theorem |2] then states criteria under 
which these sets diffract. With this result, we explain in Sec.|3]how a cut and project scheme 
naturally appears. So, we obtain the associated internal space of model sets, like the Fibonacci 
or the rhombic Penrose tilings, through the information given by the autocorrelation. Model 
sets (or cut and project sets) are the most common models for aperiodic order (quasicrystals). 

In Sec.0]we investigate lattice substitution systems. We first explore the diffraction spec- 
trum for scatterers distributed (aperiodic or even stochastic) over a lattice (Theorems|3and|5}, 
before turning to distributions obtained by substitutions. For this special case, we are again 
interested under which conditions they diffract (Theorem |6|l. As an application of a lattice 
substitution system as well as a model set with non-Euclidean internal space, we calculate, in 
Sec.|5l the diffraction spectrum of the paperfolding sequence. 
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Thereafter, we leave the area of deterministic point sets and turn to systems with disorder, 
i.e., to random tiUngs in Sec.|6l We carefully introduce the notion of a random tiling, before 
we state results about one dimensional binary random tilings in Theorem |8l and about the 
two dimensional Ising lattice gas in Theorem |9l We end this section by the example of a 
Fibonacci random tiling, where we are particularly interested in the role of the internal space 
of quasiperiodic random tilings for their diffraction spectrum (Theorems ^| and II 1> . We 
conclude our article with an outlook (Sec.0 where we indicate future directions in diffraction 
theory. 



2 Mathematical diffraction theory 

The basic object of interest is a set of scatterers in a Euclidean' space M'^, which we model 
by a translation bounded complex measure^ uo. It describes the distribution of matter in a 
mathematically adequate way. To calculate the diffraction spectrum, we need the autocorre- 
lation measure 7^ attached to uj. The diffraction spectrum is then the Fourier transform 
of the autocorrelation measure. Here, ^^{E) is the total intensity scattered into the volume 
element E, and thus describes the outcome of a diffraction experiment, compare |81. It can 
uniquely be decomposed as = {'^^)pp + (7i^)sc + (7Lj)ac ^y the Lebesgue decomposition 
theorem, see ||4J. Here, (7j^)pp is a pure point measure, which corresponds to the Bragg part 
of the diffraction spectrum, (7i^)ac is absolutely continuous and (7jj)sc singular continuous 
with respect to Lebesgue measure. The pure pointedness of the diffraction spectrum of ui 
is in question, i.e., whether the Lebesgue decomposition reduces to 7^ = {j^)pp (then, the 
diffraction spectrum consists of Bragg peaks only). 

We first fix an averaging sequence A = {-B„ | n G N} of balls of radius n around 0, so 
we begin analysing the spectrum by looking at finite pieces (balls) of the structure in question. 

We define = oj{,f), where f{x) = /(— x), and set cj„ = and u)„ = {^nY ■ Then, 
the measures 

vol(Bj ^ ^ 

are well-defined, since they are (volume averaged) convolutions of two measures with compact 
support. The autocorrelation of u exists, if (7f;"'')„gpj converges in the vague topology, 
compare (6|, and is then, by construction, a positive definite^ measure. 



'Most results also apply to more general spaces, namely tr-compact locally compact Abelian groups, with 
Lebesgue measure replaced by Haar measure, etc., see 

-Frequently, we make use of the one-to-one correspondence between measures and regular Borel measures by the 
Riesz-Markov representation theorem, where a measure is a continuous linear functional on the space of compactly 
supported continuous functions on W^, while a Borel measure is defined on the Borel sets of R**. The convolution of 
two measures ^, v is defined as * vi^f) = Jj^^ ^j^^ f(x + y) dii{x) dv(y) and is well-defined if at least one of 
them has compact support. A measure is translation bounded if, for all compact K C M'*, sup^g^^ + K) < 
< 00 for some constant which only depends on K. Here, |^| denotes the total variation measure (which is 
positive) and t + K = {t + x \ x a . See 1 2 3 4 5 6 7 and references therein for details. 

measure ^ is positive definite iff ^(/ * /) > for all complex-valued continuous functions with compact 
support, compare jSJ. 
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We say that a measure is almost periodic if the set of translates is relatively compact (i.e., 
its closure is compact). Of course, we have to fix a topology for this, and in our case we need 
the strong or product topology on the space of translation bounded complex measures (and 
not the vague topology), see |2||5). With this topology, we speak of strong almost periodicity 
to distinguish it from almost periodicity in other topologies. The key result reads: 

Theorem 1 (2||5l The measure lo is pure point dijfractive iff 

(i) 7„ exists. 

(ii) 7^ is strongly almost periodic. □ 

Note that|(i)]is merely a convention, since one can always pick an appropriate subsequence of 
the averaging sequence A for which 7^ exists. 

So, given a measure uj, we have to check these two conditions to decide whether it is 
pure point diffractive. This can be done for many relevant examples. We would also like to 
solve the homometry or inverse problem, i.e., the question which measures uj account for a 
given diffraction spectrum 7. This is a hard problem, because there is no inversion process of 
Eq. [0 In fact, rather different can have the same diffraction, see 1 9 1 . One also would like to 
understand the implications of a pure point spectrum. Here, we will not explore these last two 
questions further. 

Instead, let us specialise on the situation of a countable set S of scatterers in with 
(bounded) scattering strengths v{x), a; € 5'. It can be represented as a complex Borel measure 
in the form of a weighted Dirac comb'^ 

xes 

where is the unit point (or Dirac) measure located at x, i.e., S,^{(p) = (p{x) for continuous 
functions (p. This way, atoms are modeled by their position and scattering strengths. Convo- 
lutions with more realistic profiles are not considered here, but can easily be treated by the 
convolution theorem, see 1 3 1 . 

Denote the set of "inter-atomic distances" hy A^S — S — {x~y\x,yES}. Then we 
make the following three assumptions, see ||2|: 

(A 1) The measure co is translation bounded, i.e, there exist constants Cj^ so that 

sup ^ 1^(2^)1 < < 00 

for all compact sets K C R'^. 
(A 2) The autocorrelation coefficients 

1 



x.yeSnB^ 
x — y—z 



exist for all z G /\ (we set rj{z) = if z ^ A). Consequently (if also (A|3} holds), the 
autocorrelation measure 7^ = X^zeA '?(-^) exists. 



'^For later reference, we write to = A - 5gif v{x) = A for all a: £ S. 
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(A 3) The set A'"^^ — {z ^ A \ ri{z) 7^ 0} C Z\ is uniformly discrete, i.e., there is an ?■ > 
such that open balls of radius r centred at the points of A^'^^ are mutually disjoint. 

The support of (which is the support of rj{z)) plays a special role, in particular the 
group L generated by it: 

L = (Zi"'')^; C (2) 

The point set S is repetitive if for any compact set K C W^, {t <^W^\ SnK = {t+S)nK} 
is relatively dense^; i.e., there exists a radius R — R{K) > such that every open ball Bii{y) 
contains at least one element of t e M'' for which S* n X = {t + S){^K. If S" is repetitive, we 
have A = Zi'='% hence L = (zi)^. 

We define a translation invariant pseudo-metric^ by 



Q{s,f) = 



1) 



(3) 



If all weights v{x) are non-negative, one has g{s,t) < 1 for all s,t E M'', but, in general, 
g is bounded by V2- This pseudo-metric g defines a uniformity, both on L and K*^, see |2|. 
The induced topology is, in general, completely different from the usual Euclidean topology 
of M.'^. It is called the autocorrelation topology. 

Next we define, for e > 0, a set of e-almost periods of the autocorrelation through 

P^ = {teM.'^\git,0) <e}. (4) 

We clearly have the following inclusions for e < e': 

{periods of w} c C P^, C R'', 

and furthermore Pi = A''^^ if all weights v{x) are non-negative. Now, we are able to apply 
Theorem^to weighted Dirac combs. 

Theorem 2 |2| Let ui be a weighted Dirac comb that satisfies (A^, (A|2j and (A|3|l. Then 
exists, and the following statements are equivalent: 

(i) is relatively dense for all £ > 0. 

(ii) is norm almost periodic. 

(iii) is strongly almost periodic. 

(iv) is pure point diffractive. □ 



The norm almost periodicity in |(ii)| refers to the topology defined by the norm = 
supj.gj{d \ui\{x + K) for some fixed compact K with nonempty interior, e.g., a closed unit 
ball. The concept of norm almost periodicity is independent of the choice of K. 

This theorem applies to the diffraction from the visible lattice points, see ^ 2 , .10J . but also 
to model sets which we will consider next. 



'a set Q is relatively dense if there is a radius ij > such that every ball of radius R in contains at least one 
point of Q. 

'a pseudo-metric is a non-negative, symmetric function on R'* X that satisfies the triangle inequality. Such a 
function q is translation invariant if q(s + r,t-\- r) = g(s, t) for all r £ R"^. 
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The well-known cut and project mechanism provides examples of weighted point sets which 
satisfy the assumptions (A^, (A|3 and (A|3}- But here, we start with a countable weighted 
set of scatterers in R'^ for which the corresponding weighted Dirac comb fulfils Theorem |2] 
and is therefore pure point diffractive. This will lead us to a suitable cut and project scheme. 

The constructive picture behind this is the following: On L of Eq.|2] we have a pseudo- 
metric g, defined in Eq. |3] which in turn defines a uniformity, and then the autocorrelation 
topology on L. Inside L, the sets of Eq. |3 are the open balls of radius e around in 
the autocorrelation topology. The group L equipped with this topology admits a (Hausdorff) 
completiorP H, which is a locally compact Abelian group^. This means that there is a con- 
tinuous group homomorphism (p: L ^ H where ip{L) is dense in H. The relative denseness 
of is crucial here, because for every set of L there is an open set B{e) of H so that 
(p{P^) ~ n B{e) and B{e) has compact closure B{e) for < e < 1, giving the local 
compactness of H. 

We obtain the following cut and project scheme, see (3: 

TT 

M"* < — R'^ X H 



L < — > {{x,ip{x))\x e L} 



^ (5) 

dense 



Here, L is a lattice in M'^x H, i.e., L is a closed subgroup so that the factor group (R'' x H) /L 
is compact. The projection TTim is dense in internal space H and the projection tt into physical 
space W'' is one-to-one on L. 

The internal space H and the lattice L both arise from the group L via the autocorrelation 
topology. In spite of its abstraction, the procedure gives back the familiar Euclidean internal 
spaces of the well-known examples (e.g. Fibonacci). In the case of the rhombic Penrose 
tilings it gives the minimal internal space possible for representing them in the cut and project 
formalism: = R^ x (Z/5Z), see |2|. For other tilings, the internal space may be p-adic, 
see Section|5lfor an example. Note that the completion map is not one-to-one in general, 
its kernel is ne>o '■^^ group of statistical periods. For example, if 5 is a lattice with all 
weights equal, then S — L,ip = 0, H = and L = L,so the cut and project scheme collapses 
into triviality. In general, the internal space H ignores the periodic part of oj, for which no 
additional structure is required, and reflects only the aperiodic parts. 

A set yl C R'' is a model set for the cut and project scheme in Eq.|5j if there is a relatively 
compact setW C H with non-empty interior and a t G R'' such that 

A = t + A{W) = t + {xeL\ ip{x) e W}. 

^ H is a completion of L, if L has dense image in H and every Cauchy sequence in H has a limit in H, e.g., R is 
the completion of Q. The completion is unique up to topological isomoiphism. 

topological space is called locally compact if each point is contained in a compact neighbourhood. If this space 
is an Abehan group, we speak of a locally compact Abelian group. Examples are Euclidean and p-adic spaces. 
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Note that, in the context of model sets, the map ip is often called the *-map and denoted by 
(•)*, i.e., one writes x* = <pix)- A model set is always a Delone set, i.e., it is both relatively 
dense and uniformly discrete. A model set is called regular if the boundary of W has Haar 
measure'^ 0. Regular model sets are the most relevant model sets for the physical applications 
in the theory of quasicrystals. They also play a prominent role in the analysis of sequences 
with long-range (aperiodic) order, cf. 1,1 IJ and references therein. 
One of the cornerstones of the theory of model sets is: 

Theorem 3 l6lll2l Regular model sets are pure point dijfractive. □ 

Let us now go back to our discussion of diffraction in the context of the assumptions (AQ, 
(A|2li and (A|3|l. The pure point diffraction in Theorem|2lis intimately related to the cut and 
project scheme in Eq.|5] But it can happen that the set S itself is not a model set (e.g., as in the 
case of the visible lattice points), see |2|. So the question arises: Which pure point diffractive 
point sets are (regular) model sets? We have only partial progress on this. 

Noting that = f{Pe) has non-empty interior for all < e < 1, it follows form our 
above discussion that C yl(Wj). Since L is countable, we even get that = A{Wg) can 
be violated for at most countably many values of e, so that A'^'^^ is the union of an ascending 
sequence of model sets. Furthermore we have 

Proposition 1 fT3\ Assume (A^, (A|2j and (A|3} hold. Then A'^'^'^ is a model set. □ 

This leaves open the question of whether or not yl is a model set. Progress in this seems to 
depend on utilising the dynamical hull of A: 

X{A) = {A + t\reW\, (6) 

where closure is taken with respect to the local topology^^. Then (X(yl), R"^) is a dynamical 
system under the obvious action of K'^ on X{A). 

If yl is a repetitive regular model set then it is known 1 12 1 that X{A) is strictly ergodic, i.e., 
both minimal" and uniquely ergodic". There is a conjecture to the effect that, conversely, 
any pure point diffractive Meyer set^- of R'^ for which X{A) is strictly ergodic is in fact a 
model set. 

4 Lattice substitution systems 

An interesting class of point sets is formed by the subsets of a lattice. Even though they can be 
aperiodic or even stochastic, the underlying lattice P leaves its imprint, most notably in form 
of a periodicity of the diffraction, with the dual lattice P* as lattice of periods. 

'On every locally compact Abelian group there exists a unique (up to a multiplicative constant) translation invari- 
ant regular Borel measure. This is called the Haar measure, and is given by the Lebesgue measure on Euclidean space 

'"informally, two discrete and closed point sets are close in the local topology, if, after a small translation, these 
two coincide on a large compact region. 

"a dynamical system (X, T) is minimal in case X has no proper closed T-invariant subsets. It is uniquely ergodic 
if there is only one T-invariant Borel probability measure on X. 

Delone set S is a Meyer set iff the set of "inter-atomic distances" Z\ = S — S is also a Delone set. Every 
model set is a Meyer set. 
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In a more general formulation, let w: 7^ ^ C be any bounded function, and consider the 
weighted Dirac comb 

This includes the previous case via v = Ig, i.e., v{x) = 1 for x E S and otherwise. Then, 
if 'Y^ is any of its autocorrelations (e.g. as obtained along a suitable subsequence of averaging 
balls), we have the following result. 

Theorem 4 1141 Let F be a lattice^^ in and cu a weighted Dirac comb on F with bounded 
complex weights. Let be any of its autocorrelations, i.e., any of the limit points of the family 
{7!"'' I n £ N}. Then the following holds. 

(i) The autocorrelation 7^ can be represented as 

xer 

where O : ^ C is a bounded continuous positive definite function that interpolates 
the autocorrelation coefficients rjlx) as defined at x € F. Moreover, there exists such a 
which extends to an entire function 0: — > C with the additional growth restriction 
that there are constants C,R > and G Z such that \0{z)\ < C • (1 + \z\)^ ■ 
exp(i? I Im(z)l) for all z G C^. 

(ii) The diffraction spectrum 7^ of lu is a translation bounded positive measure that is pe- 
riodic with the dual lattice^'^ F* = {k £ R*^ | {k,x) £ "L for all x G -T} as lattice of 
periods. Furthermore, 7^ has a representation as a convolution, 

in which g is a finite positive measure supported on a fundamental domain of F* that is 
contained in the ball of radius R around the origin. □ 

This reduces the analysis of the spectral type of 7^ to that of g, which has compact support. 
An interesting application concerns lattice subsets and their complements, which leads to 
the following result. 

Theorem 5 II14I Let F be a lattice^^ in M'', and let S CZ F be a subset with existing (natural) 
autocorrelation coefficients Tjg{z) — dens (5 H (z + S)). Then the following holds. 

(i) The autocorrelation coefficients rjga{z) of the complement set S'^ ~ F \ S also exist. 
They are rjgc (z) = for all z ^ F and otherwise, for z £ F, satisfy the relation 

ris<:{z) — dens(S"^) = Vsi^) ~ dens(S') 

" These results can be generalised to lattice subsets in locally compact Abelian groups whose topology has a 
countable base, see 1 14 for details. 

''*The Euclidean scalar product is denoted by (■, ■). 
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(ii) If, in addition, dens(S') = dens(-r) /2, the sets S and S'^ are homometric. 

(iii) The diffraction spectra of the sets S and S'^ are related by 

■jg, = 75 + (dens(S"') - dens(S')) • dens(r) 6^, . 
In particular, 75c — 75 if dens(S''^) = dens(S'). 

(iv) The diffraction measure 7^0 is pure point iff^g is pure point. 



□ 



As an immediate consequence of part (i), one can check that the two pseudometrics defined 
by S and S""^ via Eq.|3lare scalar multiples of one another, hence define the same uniformity 
(as long as dens(S') • dens(5'') > 0). 

Of considerable interest in this context are lattice subsets which are obtained by lattice 
substitution systems via Delone multisets^^. They are the natural generalisation of one dimen- 
sional substitution rules with constant length, and it was recently possible L15.. .16. .17 I to find a 
complete generalisation of Dekking's coincidence criterion (see 1181 ^ to this general situation. 
The result, stated below, is a circle of equivalences which directly puts pure pointedness and 
model sets on the same footing. Although one of the equivalent criteria is modular coinci- 
dence, we refer the reader to 1171 for the rather technical definition, just pointing out that its 
primary virtue is that it is testable by a straightforward algorithm. 

A matrix function system (MPS) on a lattice C M"* is given by an to x m-matrix 
<P = (^y )' where each <P^j is a finite set (possibly empty) of mappings F ^ F. Here, the 
mappings of <P are affine linear mappings, where the linear part has the form x 1^ Q x and 
is the same for all maps. We call Q the inflation factor. It is required to have all eigenvalues 
exceeding 1 in absolute value. Any MPS induces a mapping or substitution on P{F)™, 
where P{F) denotes the set of subsets of F: 



' U U f{u^) \ 



, U U f{u^) 



(7) 



We say that U — {Ui, . . . , is a fixed point of ^ if [/ = <1>U. Purthermore, we call 
{U, <P) a lattice substitution system on if is an MPS on F, U is a fixed point of the 
U/s are pairwise disjoint and all the unions in Eq.0are disjoint. A lattice substitution system 
is primitive if its corresponding substitution matrix M — {4/^'l'ij)i<:i j<,„ is primitive, i.e., if 
there is a fc € N such that Af'^ has positive entries only. 

Let U he a lattice substitution system for the MPS We say that a (7-cluster C = 
U n Bj^{s) — (t/j n Bj^{s))^^j^, defined by intersecting all the components of U with a 

'^A multiset in M"^ is a subset U-^ X ■■ ■ X C M'* X ■ ■ ■ X R'* (m copies), where C K''. We also write 
U= ([/j^, . . . , C/^)* = (C^i)i<„- We say that (!7j)^^^ is a Z)e/one muW.vef in R** if eacli !7. is Delone and the 
union UIILi C K"* is also Delone. It is convenient to think of a multiset as a set with types of colours (types of 
atoms), i being the colour of points in U^. 
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ball of radius R around a point s E M.'^, is legal if it lies in <^"(tt) for some point u of U. 
Furthermore, we define the symmetric difference of two multisets as symmetric difference of 
their corresponding components, i.e., UAV = ([/,AFJ,<,„ = (((7, \ V,) U {V, \ t/,)),<„,. 
We also use the notation S + T = {x + y \ x e S, y e T} and U + S ^ + S')i<„ = 
({a; + y \ X e U^, y e S})^^„^ for sets S, T and a multiset U. 

Theorem 6 1171 Let U be a primitive lattice substitution system with expansive map Qfor the 
lattice r — Ui<m in^'^ cind suppose that every U -cluster is legal. Let F' = ri + - ■ ■+r^, 
where = {A^ — A^)-^. The following assertions are equivalent. 

(i) U has pure point diffraction spectrum (meaning that each has this property); 

(ii) U has pure point dynamical spectrum (meaning that each has this property); 

(iii) dens([/A(Q"a + U)) = (0),<„/or all a G F' ; 

(iv) A modular coincidence relative to Q^^ F' occurs in 'P'^^ for some M; 

(v) Each Uj^ is a regular model set, i < m, for the cut and project scheme 

^ — Rrf X r — > T 
u 

F < — f 

□ 

Here, the pure point dynamical spectrum is defined as follows: For a dynamical system 
{X, T) (see Eq.|6|l that has a (unique) invariant probability measure /i associated to it (which is 
the case for dynamical systems which arise from primitive substitutions), we have the Hilbert 
space L'^{X,^) and the unitary operator B: L^{X,^) L'^{X,^), f ^ f ° T. If the 
eigenfunctions of B span L'^{X, fi), then we have a pure point dynamical spectrum. Note that 
the equivalence of|(i)land |(ii)| can also be established in a more general setting, see |,19.| . 

The model set in this theorem is with respect to a very particular cut and project scheme. 
The key point is the internal group F which we now briefly explain. F is the Q-adic comple- 
tion (in terms of a profinite group, see L20u21J ) 

T = (r)n = lim f/q'^f' = iim(r/r' ^ f/qf' ^ . . . ^ f/q'^f' ^ ...) 

of F, supplied with the usual topology of a profinite group (which makes it compact). We 
note fliat F embeds naturaUy into T. Then, F is the group {{t, t) e R"^ x T 1 1 e T}. 

This gives a satisfactory approach to those systems which are model sets, including an 
algorithm to test it. The latter, however, is of limited value to disprove the model set property, 
unless one can see that no coincidence can ever occur For this situation, Frettloh has recently 
proved several sufficient criteria. Although they are not exhaustive, they are easy to check and 
seem to cover many cases of relevance, see L22J for details. 
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5 Paperfolding sequence as model set 

To demonstrate the usefulness of the general setting, let us consider an explicit example with 
a p-adic internal space. 

The so-called regular paperfolding sequence'^ starts as 

11011001110010011101100011001001 ... (8) 

and can be obtained by folding a sheet of paper repeatedly to the left, see ll23l : 



X 



The sequence in Eq. |8l is obtained by unfolding the (infinite) stack and encoding a left 
(right) bend by 1 (0). 

An alternative description employs two steps: The first determines the unique one-sided 
fixed point of the primitive 4-letter substitution of constant length'^: 

a — > a6 
b ^ cb 
c ^ ad 
d — > cd, 

which starts as abcbadcbabcdadcb .... The following second step maps a and 6 to 1 and c and 
d to 0, giving the sequence in Eq.|8] 

Here, we change the point of view slightly, in that we consider two-sided (or bi-infinite) 
fixed points of a, of which there are precisely two'^: 

b\a —>■ cb\ab — > adcb\abcb —>■ ... — > Wi 
d\a cd\ab adcd\abcb ... — > W2 

where | denotes the seamline. Note that Wi and differ only in the first position to the left 
of the seamline, otherwise they are equal. 

Let us represent the letters with intervals of equal length 1, and points on their left end, of 
type a, b, c and d. Let denote the a-points, etc. Then, the substitution together with the 



More generally, we can define a paperfolding sequence recursively: the sequence {01,02, . . . } is called a 
paperfolding sequence iff ai = —03 = as = —07 = . . . , and the remaining sequence {02, 04, . . . } is a paper 
folding sequence. We obtain the regular paperfolding sequence for a„ £ {1,-1} and writing for —1. 

"This is a primitive substitution of height 1 which has a coincidence (after two substitutions). Therefore it is pure 
point diffractive by a criterion of Dekking, see 1 18|. 

'**The dynamical zeta function can be calculated with the method of Anderson and Putnam 1 24 1, and gives 



l-2z 



From here, one sees that the number of fixed points of ct'" (m > 1) is 2™, which shows that w-^^ and the only 

solutions of (T(-!ii) = w with w a bi-infinite word. 
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fixed point equation a{w) — w lead to the following system of equations: 

2i7„ U 2i7^ 

= (2^?a + l) U (212^ + 1) 

/2rf = (2^2^ + 1) U (212^ + 1) 

where also il^ U i7j U U = Z by construction. With this, one quickly checks that the 
first and the third equation lead to the unique solution 

I2„=4Z, i7^ = 4Z + 2 

which reduces the other equations to 

= (8Z+1) U (2 12^ + 1) 
^ (8Z + 5) U (2i7^ + l). 

Since and are subsets of Z, the general solution is 

^ [j 2"+2 Z + 2™ - 1 

m>l 

= U 2™+2Z + 3-2'"- 1, 

m>l 

where the only remaining freedom consists in adding the singleton set { — 1} to either of them. 
This reflects the difference between the two fixed points, Wi ({ — 1} goes to J7j,) and W2 ({^1} 
goes to n^). 

If one now follows the construction of a canonical cut and project scheme as derived in j2|, 
one finds that the autocorrelation topology is the 2-adic topology, and this completes Z (the 
set of differences) to Zj, the compact group of 2-adic integers. So we have the cut and project 
scheme 

R < — M X Z2 — > Z2 
U 
L 

where L = {(m, m) | m e Z} is a lattice. Our points are now model sets in this scheme. 
Defining the windows (as subsets of Z2) 

one finds n = {— 1} and thus 

n,^A{w,), n,^A{w,), n,^A{wj, n^ = A{wA{-i}) 



for Wi, while { — 1} moves from to for u>2. 
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As a regular model set, the paperfolding sequence is pure point diffractive'''. More pre- 
cisely, if Lu = A ■ + ^ ' '^12 + ^ ■ ^r2 + ^ ■ '^f2 ' '■^^ diffraction measure reads 



A+B+C+D 



E 



A- B+C -D 



A~C 



E 

■r>3 



Finally, let us consider the binary reduction. One gets 

n^un^ = u 2"+2 z + 2™ - 1 

U 



m>0 

U 2'"+2 Z + 3 • 2" - 1 

m>0 



(9) 



plus { — 1} added to one of them. Clearly, these are again regular 2-adic model sets, and thus 
also pure point diffractive. To summarise: 

Theorem 7 The quaternary regular paperfolding sequences Wi and are regular 2-adic 
model sets, with pure point diffraction spectrum as given in Eq. |9] Also, the binary reduction 
is a regular 2-adic model set, hence also pure point diffractive. □ 

This structure is then inherited by the entire Ll-class^" of the paperfolding sequence. The 
members can be obtained via different folding sequences^', see |23 25 1 and references therein 
for details. Further examples along similar lines can be found in L2,. .26..27J . 



6 Systems with disorder 

Here, we consider diffraction properties of stochastic point sets. Simple, well-understood 
examples are Bernoulli subsets of lattices or model sets 1 141 1281 . and certain lattice gases, 
which can be analysed using elementary methods from stochastics. This approach has recently 
been generalised considerably | 29| to cover stochastic selections from rather general Delone 
sets. The results prove the folklore claim that uncorrected random removal of scatterers 
has two effects, namely reducing the overall intensity of the diffraction of the fully occupied 
set, without changing the relative intensities, and adding a white noise type constant diffuse 
background. The influence of disorder due to thermal fluctuations is discussed in 1301 . 

"it is an example of a limit periodic system. 

^"Two structures and ylj are locally indistinguishable (or locally isomorphic or LI) if each patch of A-^ (es- 
sentially, the intersection of A-^ with a compact set) is, up to translation, also a patch of /Ij and vice versa. The 
corresponding equivalence class is called Ll-class. 

^'We get the paperfolding in Eq.|8]by setting 1 = = = = . . . = = . . . (see Footnote ll6l . because 
we only fold in one direction. Using a different folding sequence, which corresponds to the aj,, not being equal, 
we get a different paperfolding sequence, but this one is LI to the one in Eq.lUand vice versa and therefore in the 
same Ll-class. So, all such paperfolding sequences have the same diffraction spectrum, namely Eq.|9]in the binary 
reduction, i.e., A = B and C = D. 
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A prominent class of stochastic point sets are random tilings 1311 1321 l33l . Diffraction 
properties of these tilings are understood for systems without interaction, for one-dimensional 
Markov systems, and for product tilings 1 91 1341 l35l . Systems with interaction are generally 
difficult to analyse. Here, only few results are available for certain exactly solvable models 
from statistical mechanics with crystallographic symmetries, whose autocorrelation can be 
computed explicitly i9l 1341 1351 1361 . Symmetries of a stochastic point set are understood to be 
symmetries on average. For a detailed discussion of this concept, see 1321 1341 . 

Most examples of random tilings with quasicrystallographic symmetries are obtained from 
ideal quasicrystallographic tilings by relaxing the allowed local configurations. Since, as for 
ideal tilings, random tiling coordinates may be lifted into internal space via the *-map, the 
random tiling ensemble possesses a so-called height representation which, in a way, can be 
understood as a description on the basis of a deviation from a model set. At present, there are 
no rigorous results concerning diffraction properties of random tilings with height representa- 
tion in dimension d> 2. Henley 1 3 1 1 argues, using elasticity theory for the free energy of such 
an ensemble, that the discrete part of the diffraction spectrum only consists in the trivial Bragg 
peak at the origin in dimension d < 2, since the width of the distribution of scatterer positions 
in internal space diverges with the system size. In dimensions d > 2, the distribution width 
converges with the system size, implying a non-trivial discrete part in the diffraction spectrum. 
Even less is known about the nature of the continuous part in d > 1, though absolute conti- 
nuity is expected. For a numerical investigation of diffraction properties of the randomized 
Ammann-Beenker tiling, see [35 1. In what follows, we will focus on stochastic disorder in 
random tilings, mainly in the one-dimensional case, because the understanding of the higher 
dimensional situation is still rather incomplete. 

The diffraction of ID random tilings^^ has been investigated previously I34il . ID binary 
random tilings have a non-trivial pure point part iff they have a rational interval length ratio 
a = u/v. 

Theorem 8 11341 The natural density of A exists with probabilistic certainty and is given by 
d — [pu + qv)~^. If uj — 5y^ — X^xeA denotes the corresponding stochastic Dirac comb, 
the autocorrelation of lu also exists with probabilistic certainty and is a positive definite 
pure point measure. The diffraction spectrum consists, with probabilistic certainty, of a pure 
point (Bragg) part and an absolutely continuous part, so 7^ = (7^)pp + (7(j)ac- 
If a = u/v, the pure point part is 



where, if a € Q, we set a — a/b with coprime a, 6 G Z and define ^ — uj a — v/b. The 
absolutely continuous part can be represented by the continuous function 



^^A (ID binary) random tiling [31 32 is a covering of the real line with two intervals of fixed lengths u and v 
without gaps or overlaps. Associated with each random tiling is the set A of left endpoint positions of its intervals. 
We call two random tilings equivalent if they are equal up to a translation. For each equivalence class, we choose a 
representative with S yl. The random tiling ensemble is the set of all non-equivalent random tilings. 




9{k) 



d ■ pq sin^(7rfc (m — u)) 



p sin^ (nk u) + q sin^ (nk v) — pq siv? {irk (u — v)) 
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Figure 1: Absolutely continuous background of a Fibonacci random tiling. 



which is well defined for k {u 
If a is irrational, this is g{k) = 



- v) ^ Z. It has a smooth continuation to the excluded points. 
Q for k{u — w) e Z with k ^ and 



5(0) 

For a — a/b Cz 
kv (f. Zj, and 



d ■ pq{u ~ w)^ 



9ik) 



Q as above, it is g{k) 
(pa + qb)^ 



{pu + qvY 

for fc (w — w) e Z, but ku ^ Z (or, equivalently. 



for the case that also ku 1,. 



□ 



The most prominent ID random tiling is the Fibonacci random tiling^^. According to be 
above theorem, its diffraction spectrum consists with probabilistic certainty of a trivial Bragg 
peak at the origin and of an absolutely continuous background, see Fig. The absolutely 
continuous background shows localised, bell-shaped needles of increasing height at sequences 
of points scaling with the golden ratio r. This is reminiscent of the perfect Fibonacci tiling. 

In dimensions d > 2, properties of the autocorrelation are known only for certain simple 
systems of statistical mechanics with crystallographic symmetries which can be interpreted in 
terms of dimer systems, see 1341 1361 l37l . This includes the domino and the lozenge tiling, the 
Ising lattice gas, and others. Here, the asymptotic behaviour of the autocorrelation coefficients 
can be computed explicitly, leading to proofs of existence of an absolutely continuous part in 
addition to a pure point part. 



^'A Fibonacci random tiling is a random tiling with interval lengths u = t = {1 + \/5)/2 and v = 1, with 
occupation probabilities p = 1/t and g = l — p=l/T^of the intervals (almost surely). Each interval endpoint of 
a representative of a Fibonacci random tiling belongs to the module Z[r] = {mr + n | m, n £ Z}. Every (ideal) 
Fibonacci tiling also appears as a Fibonacci random tiling. 
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Theorem 9 1341 Away from the critical point, the diffraction spectrum of the Ising lattice gas 
almost surely exists, is 1? -periodic and consists of a pure point and an absolutely continuous 
part with continuous density. The pure point part reads 



where the density p is the ensemble average of the number of scatterers per unit volume. At 
the critical point, the diffuse scattering diverges when approaching the lattice positions of the 
Bragg peaks. 

We now analyse the effect of a cut and project setup on diffraction properties. For ID 
random tilings, an embedding into is given as follows. Each position a; G yl of an element 
of the random tiling ensemble may be written in the form x = mu + nv, where m, n G Z. 
If a = u/v ^ Q, the numbers m, n are uniquely determined. If a e Q, uniqueness is 
achieved by parametrising a; = Obym = ri = and incrementing (decrementing) m by 
addition of a w-interval to the right (left), and likewise with n for t;-intervals. Identifying the 
unique coordinates m, n with points (m, n) S 1?, we map a random tiling to a bi-infinite, 
directed walk on the edges of the square lattice. We call two walks equivalent if they are 
equal up to a translation. For each equivalence class of walks, we may choose a representative 
which passes through the origin in 1?. This establishes a one-to-one correspondence between 
non-equivalent random tilings and non-equivalent bi-infinite, directed walks. 

Let us restrict to Fibonacci random tilings. Recall that the (ideal) Fibonacci tilings may 
be obtained within the cut and project setup by (orthogonal) projection of lattice points of a 
scaled copy of 1? confined to a suitable strip onto the subspace with irrational slope 1/t. This 
way, we obtain a one-to-one correspondence between (random) tiling coordinates in direct 
space and in internal space via the *-map, given by (tot + nY — mr' + n, where r' = —1/t 
is the algebraic conjugate of r. The value x* of a tiling coordinate x is also called its height, 
and the above collection of direct and internal space together with the canonical projections is 
also called height representation of the Fibonacci random tiling ensemble. 

In the following, we will consider the distribution of scatterer positions of Fibonacci ran- 
dom tilings in internal space. We restrict ourselves to patches of Fibonacci random tilings of 
M consecutive intervals on the positive half-axis, starting at a; = 0. Following |351, we con- 
sider the occupation probability for the position of the rightmost interval. Since the random 
tiUng patch is a Bernoulli system, the probability of the position being x = rriT + [M — rn), 
or equivalently x* = mr' + (Af — m) in internal space, is given by 



According to the theorem of de Moivre-Laplace, the binomial distribution (with mean fi — 
Mp and variance cr^ = Mpq) may be approximated by the Gaussian distribution for large 
fixed M, yielding 



1. T>Tc: {%)pp = 
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Note that in the limit M oo, the admissible positions x* lie dense in internal space. We 
fixed the normalisation such that the integral of the density over internal space equals unity. 
To obtain the distribution of all interval positions of random tiling patches of length N, we 
sum over all endpoint positions of patches with n < N intervals and normalise^'*, 

1 ^ 1 /"^ 

p^^*^ = ^ E - 1^ / p(''' (10) 

In the second equation, we approximated the sum to leading order in N using the Euler- 
Maclaurin summation formula. This leads to the following result. 

Theorem 10 The distribution p{x*) of scatterer positions of Fibonacci random tiling patches 
in internal space is, to leading order in the patch size N, given by 

where erfc(x) = e~* dt denotes the complementary error function. □ 

The function f{z) is shown in Fig.|2] Note that the width of the distribution grows with 
the system size as ^/N, but the distribution itself is not a Gaussian, as usually believed^^. 

This result was derived for patches of intervals, starting from the origin in positive 
direction. Since the result is symmetric in x*, it is also valid for tiling patches with N intervals, 
starting from the origin in negative direction, hence also for tiling patches with 2N intervals 
with N intervals on the negative and on the positive axis each. 

more natural normalisation may be the density of points instead of unity, see below. 
^^For a comparison with numerical simulations of the Fibonacci random tiling ensemble and for the example of 
the two-dimensional Ammaim-Beenker random tiling, see ^35t . 
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Within the cut and project scheme, the Dirac comb of a model set may be characterised as 
a sum of point scatterers over the projected lattice points 



weighted by the characteristic function of the (compact) window in internal space. For diffrac- 
tion of quasicrystallographic random tilings, the common approach I31J is to investigate prop- 
erties of an averaged structure given by the Dirac comb weighted by the occurrence probabil- 
ity of each scatterer within the random tiling ensemble. The averaged distribution in internal 
space need not exist for the infinite tiling, as we showed in the previous section. It is generally 
believed 1311 that it exists in dimensions d > 2, the case d = 2 being marginal with a log- 
arithmic growth of the distribution width with the system size, leading to the statement that 
there is no non-trivial discrete component in the diffraction spectrum for d <2 (for aperiodic 
systems). 

As argued in Eq.^| the averaged structure may be written in the form 



where the support of ip is generally the whole internal space. Note that this object is generally 
ill-defined, since summation is over a dense set. In averaging, one loses information about 
correlations between scatterers, so that the analysis will at most yield information about the 
discrete part of the diffraction spectrum, but not about the continuous part, see the theorem 
below. We investigate under which conditions Eq. [^defines a tempered distribution. To this 
end, we assume that Lp vanishes sufficiently rapidly at infinity, which includes the Gaussian 
(and, in a certain sense, characteristic functions). Following |35|, we consider the special 
situation of a Euclidean internal space H = K™ and assume that the canonical projections tt 
and TTint are both dense and one-to-one. We denote by yo\{FD) the volume of a fundamental 
domain of the lattice L eR"^ x H w.r.t. the product measure of the Lebesgue measures on R'' 
and on H. We denote the dual lattice of i by (i)* = {x e M.''- x H\{x,y) eZ for all y e L}, 
and its projection by L* = tt{{L)*). 

Theorem 11 1381 Assume ip : M™ C continuous and limj^^oo ip(y) = for 

some a > 0. Then, the weighted Dirac comb in Ea. \ll\ is a translation bounded measure. It 
has the unique autocorrelation 



being a translation bounded, positive definite pure point measure. Its Fourier transform is a 
positive pure point measure. If it and TTint ore orthogonal projections, it is explicitly given by 




(11) 





where ip denotes the Fourier transform of ip. 



□ 
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A natural choice for the normalisation of the function ip arises form the observation that for 
Dirac combs satisfying the assumptions of the above theorem the density of points p exists. 



The above theorem sheds light onto the example of the Fibonacci random tiling. Here, 
the internal distribution may be described by a sequence of distributions of increasing width 
but constant mass. The limit of the corresponding sequence of measures will have a trivial 
discrete part and a continuous part, whose form may be compared to the above results. 



For our discussion of pure pointedness, we made substantial use of the uniform discreteness of 
^ess rpjjjg ^^jj certainly be relaxed, as Theorem[2shows, but things become considerably more 
involved beyond this "barrier". This is also intimately related to stepping into the territory of 
mixed spectra, which seems particularly timely. 

Pure pointedness of the diffraction is equivalent to strong almost periodicity of the au- 
tocorrelation. More generally, one can show that any autocorrelation possesses a unique 
decomposition into a strongly almost periodic part and a weakly almost periodic part with 
zero volume average, compare |5|. So, we get 

where "weak" refers to the weak topology in relation to the strong (product) topology. 

The Fourier transform of (7j^)sap is a pure point measure, while that of (7^)o-wap is con- 
tinuous 1 5 1, so that one has full control of this question on the level of the autocorrelation. 
Also, important issues of the diffraction of random tilings can be formulated and understood 
in this context, but most results are folklore, and still need to be proved. Furthermore, there 
is no such decomposition known that would allow a distinction of absolute versus singular 
continuity. It is highly desirable to improve this situation in the future. 

Finally, even if all the spectral questions were settled, the big remaining question is how 
to characterise the homometry class, i.e., the class of measures with a given autocorrelation 
(and hence diffraction). This is part of the inverse problem, where results are rare at present. 
So, the question of the title should now be replaced by another one: 
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"Which distributions of matter are homometric?' 
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